Fidelity and quantum chaos in the mesoscopic device for the Josephson flux qubit 
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We show that the three-j unction SQUID device designed for the Josephson flux qubit can be 
used to study the dynamics of quantum chaos when operated at high energies. We determine the 
parameter region where the system is classically chaotic. We calculate numerically the fidelity or 
Loschmidt echo (LE) in the quantum dynamics under perturbations in the magnetic field and in the 
critical currents, and study different regimes of the LE. We discuss how the LE could be observed 
experimentally considering both the preparation of the initial state and the measurement procedure. 
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Ultrasmall Josephson devices have been used as tools 
for studying quantum phenomena at the macroscopic 
level since the 1980s Macroscopic quantum tunnel- 
ing Q and macroscopic quantum coherence of flux Q and 
persistent currents (J| have been observed experimentally. 
More recently, mesoscopic Josephson devices have been 
used for the design of qubits for quantum computation 
[S B 0> 0]- The progress made in this case allows to 
have nowadays Josephson circuits with small dissipation 
and large decoherence times, giving place to a coher- 
ent manipulation of the system 0, Q. This could also 
make possible the use of Josephson devices for the study 
of the quantum dynamics of chaotic systems, a subject 
where there has been a great interest in the last years 

SESEIEiElEl. 

The stability of the quantum dynamics of a system 
against perturbations can be quantified by the fidelity 
or Loschmidt echo (LE) [l0| . The LE is the overlap be- 
tween two states that evolve from the same initial wave 
function under two slightly different hamiltonians, 

F{t) = \f(t)\ 2 = \^ \e lH ^ R e- m ^ h \^ )\ 2 , (1) 

with Hq, H s the unperturbed and perturbed hamiltoni- 
ans, respectively. In classically chaotic systems, above 
a perturbative regime where the LE has a Gaussian de- 
cay for short times, the LE shows an exponential decay 
for large t, F(t) cx [Iq|, EI El, El El ■ For weak 

perturbation strength e the decay constant 7 depends as 
7(e) oc e 2 , and is obtained from the Fermi golden rule 
(the FGR regime) [ll| , while for strong perturbations 7 
becomes independent of e and saturates at the classical 
Lyapunov exponent A (the Lyapunov regime) • There 
have been some recent experimental measurements of the 

le [H,EiE3- 

Here we will show that the device for the 
Josephson flux qubit (DJFQ) studied in @, 0, @] is also 
a promising system for the experimental observation of 
the LE. 

The DJFQ consists of three Josephson junctions in a 
superconducting ring Q that encloses a magnetic flux 
<J) = /$0j with $0 = h/2e. Two of the junctions have the 
same coupling energy Ej and capacitance C, while the 
third junction has couplings aEj and aC, respectively 



(0.5 < a < 1). Typically the circuit inductance can be 
neglected and the phase difference of the third junction 
is: (fi3 = 2tt f + ip 1 — if 2, leading to the Hamiltonian [6( 



U = -P M P + EjV(if) 



where tp = (ipi, </? 2 ), P — M ■ df/dt, and 

M = [ M 2 fl + a + 7 -a \_ & 
' 2icJ \ -a 1 + a + 7 J v 2 Ej 



(2) 
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Here 77 = y/8E c /Ej with E c = e 2 /2C, we include in M 
the on-site gate capacitance C g — 7C, and 

V(<f) = 2 + a~cosifi~cosf2 — ctcos(2irf + ifi—if2) (3) 



In the quantum regime, P — — iKS7 v — — 2^(7^ 
the time-dependent Schrodinger equation is 
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where we normalized time by t c = H/tjEj, energy by 
Ej and momentum by h/rj. We see in Eq.(4) that the 
parameter 77 plays the role of an effective H. For quan- 
tum computation implementations @, 0, Q the DJFQ 
is operated at magnetic fields near the half-flux quan- 
tum (/ = 1/2 + 5f). In this case the two lowest energy 
eigenstates are symmetric and antisymmetric superpo- 
sitions of two states corresponding to macroscopic per- 
sistent currents of opposite sign. These two eigenstates 
are energetically separated from the others (for small Sf ) 
and therefore the DJFQ has been used as a qubit @, 0, [1] . 
As we will discuss here, the higher energy states of the 
DJFQ show quantum manifestations of classical chaos. 

It has been found in [3, E3] that superconducting 
loops with three Josephson junctions and on-site capac- 
itances (7 = 00) can be chaotic. Most Josephson cir- 
cuits, like the DJFQ, have small on-site gate capacitances 
(7 ~ 10~ 2 ). Here we analyze the dynamics of the DJFQ 
considering the realistic case with 7 = 0.02 [6(. We first 
obtain the classical dynamical evolution integrating the 
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FIG. 1: (a) Chaotic volume v c h and average Lyapunov ex- 
ponent A versus energy E for a = 0.8 and / = 1/2. Inset: 
surface of section ip xy = {fi-(p 2 )/V2, P xy = (Pi-P 2 )/v / 2 at 
£ = 1.6Ej. (b) Energy boundaries for the regimes of regular 
orbits, soft chaos and hard chaos as a function of 5f = f — 1/2 
for a = 0.8. The continuous line corresponds to the potential 
energy minimum, Emin- Inset: Energy boundaries for the 
different regimes as a function of a for / = 1/2. 



Hamilton equations that correspond to Eq.(2) with a sec- 
ond order Verlet algorithm. For different values of the 
parameter a and magnetic field / we compute the max- 
imum Lyapunov exponent A for each orbit at different 
energies E. We estimate the chaotic volume v c h(E), as 
the probability of having a chaotic orbit (i.e. A > 0) 
for a given E, using 10 3 initial conditions randomly cho- 
sen with uniform probability within the available phase 
space. Also the average Lyapunov exponent, X(E), of the 
chaotic orbits is obtained. In Fig. 1(a) we show v c h(E) 
and X(E) for / = 0.5 and a — 0.8. Above the minimum 
energy of the potential, E m i n , we find: (i) regular orbits 
for E min < E < E ch (v ch = 0), (ii) soft chaos (i.e., co- 
existence of regular and chaotic orbits, < v c h < 1) for 
E c h < E < E) lc , a Poincare section for this case is shown 
in the inset of Fig. 1(a), and (iii) hard chaos (all orbits are 
chaotic, w c h = 1) for E > E c h- The average Lyapunov 
exponent is A > above E c h- In FigQJb) we show the en- 
ergy boundaries of the different regimes (E m i n , E c h, Ehc) 
as a function of 6f = f — 1/2 for a = 0.8. We can see 
that for / = 1/2 the onsets of soft and hard chaos, E c h, 
Ehc, are closer to £ ml „ in comparison with other values 
of /. In the inset of FigQJb) we also show the boundaries 
of the dynamic regimes as a function of a for / = 1/2. 
The hamiltonian dynamics of Eq.(2) are a good approx- 
imation of the problem for energies E < 2A, with A the 
superconducting gap. The Al/AlO^/Al junctions of @ 



have 2 A w 3.7Ej. Since we find Ef lc — 1.75Ej, there is a 
wide energy range where the hard chaos regime is experi- 
mentally accessible. Furthermore, we find that for realis- 
tic experimental parameters (rj = 0.1 — 0.5, a = 0.7—0.8), 
the third quantum energy level can be above the classi- 
cal onset of chaos, E c h, for / = 1/2. In Ref. for 
the case with 7 = 00, an analysis of the statistics of the 
energy spectra of the quantum hamiltonian shows that it 
belongs to the gaussian orthogonal ensemble in the hard 
chaos regime. We find a similar result in our case, with 
7 = 0.02 and / = 1/2, for E > E hc 




FIG. 2: The average Loschmidt echo, F(t), as a function of 
time for different values of the perturbation strength of (a) 8a 
and (b) Sf. The energy is E ~ 3Ej, deep in the hard chaos 

region, and h e s = V — \J = 0.17. Time is measured in 

units of tj = h/Ej = i]tc- Insets: time evolution of a 2 F — 
2 

F 2 (t) — F(t) . The dashed lines represent an exponential 
decay given by the Lyapunov exponent (A = 0.182 ± 0.008) 
with the t = offset shifted for clarity. 

We now calculate the quantum dynamics of the D JFQ 
integrating numerically Eq. (4) with a fourth order split- 
operator algorithm as in [21j, with a discretization grid 
of Atp = 2vr/128 and At = 0.015r?t c . We use 27r-periodic 
boundary conditions on tp = (tpi,<p2). To compute the 
LE of Eq.(l), the simulations are started from minimum- 
uncertainty 27r-periodical wave packets [22[ given by 

1*0) = ^k f o {0) = C ■ e^oiv-Vo) . e -B(<p- V - )/2a 2 ^ 

with B(ip) — 2 — cos ipi — cos ip 2 , [B(<p) w \ip\ 2 /2 for small 
\ip\], Kq = (fci, k 2 ) with k%, k 2 integers, and a is the width 
of the wave packet. The LE is usually computed as the 
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average F(t) over different |\l/o) (see [12], LLiJ). Here we 
average over 15 initial conditions with different Ko,ipo 



corresponding to the same classical energy. 

T 2 



We choose 



a" = 0.3177, which corresponds to a spectral width of 
AE w 0.3Ej. We consider E = 3Ej > E hc , for which 
the classical phase space is filled by a connected region 
of chaos with Lyapunov exponent A = 0.182 ± 0.008. 
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FIG. 3: (a) Decay rate of the LE, 7, as a function of the 
perturbation strength Sf for different values of 77, with E ~ 
3Ej. The classical Lyapunov exponent is indicated in the 
dashed line. Inset: crossover perturbation 8f c as a function of 
77 _1 . (FGR: Fermi golden rule regime, LR: Lyapunov regime), 
(b) Same as (a) for perturbations 8a. 



We evaluate the LE of the DJFQ against perturba- 
tions in the parameters a and /. In FigOja) we show 
the time dependence of F(t) for different perturbations 
in the parameter a = ao + 6a for cvo = 0.8, while in 
Fig. [DJb) we show F(t) for different perturbations in 
/ = 1/2 + 6f. We can see that in both cases F(t) decays 
with time, and that the decay rate tends to increase when 
increasing the perturbation. Above a crossover value of 
the perturbation (6a c or 6f c ) we find that the curves of 
F(t) tend to overlap. The overall decay of the LE is 
has been described with the form F(t) ~ Ae~ 7t + F^, 
[13]. For large perturbations the decay rate 7 satu- 
rates at a value close to the Lyapunov exponent A of 
the classical dynamics. This can be seen in Fig. 2 where 
the dashed lines show the slope of a decay rate with 
the Lyapunov exponent for comparison. The constant 
value Faa is proportional to the inverse of the fraction of 
the volume of the Hilbert space spanned by the initial 
wave function |<3/ ) [l3|- in our case this corresponds to 
cx a 2 /(2n) 2 = 0.3l7y/(27r) 2 « 0.0013 for 7/ = 0.17, 
which is close to the results of Fig. H(a),(b). An analysis 



of the variance of the fidelity, a 2 F (t) = F(t) 2 — F(t) is 
also important [3]. This is shown in the insets of Fig. [5] 
We find that for increasing perturbations <j 2 f saturates to 
a decay given by a F (t) ~ e~ 2Xt , as discussed in 

We obtain the decay rate 7 fitting the exponential form 
Ae -7 * + Foo for F(t) for times above the initial gaussian 
decay (we have chosen t > 15tj in this case). Fig. [3£a) 
shows the obtained 7 as a function of the perturbation 
6f for different values of 77. For small perturbations we 
obtain a quadratic law dependence of the decay rate with 
the perturbation strength, 7 cx (6f) 2 , which corresponds 
to the Fermi Golden Rule (FGR) regime [n]]. F° r large 
perturbations the obtained values of 7 have a large error, 
which is of the size of the oscillations seen in the data in 
Figl^a) for large 6f. However, when comparing different 
cases of 77, we see that for large perturbations the values of 
7 fall close to the Lyapunov exponent (shown as a dashed 
line). We obtain an estimate of the crossover value 6f c 
where 7 saturates to A. In the inset of Fig. [3{a) we see 
the dependence of 8f c with 77. We find that 6f c decreases 
for decreasing 77 since quantum fluctuations become less 
important, and therefore the classical Lyapunov decay is 
reached more easily. We find a similar behavior for the 
dependence of 7 with 6a and 77 [shown in FigfSJb)] . 
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FIG. 4: Schematic setup for the observation of the Loschmidt 
echo in a DJFQ. (a) Circuit for preparation of the initial wave 
packet. V s is a voltage source and the SQUID connects a 
node of the DJFQ with a superconducting reservoir. Only 
one voltage source and one SQUID are drawn, for simplicity, 
(b) Layout for the measurement of the Loschmidt echo, ch- 
DJFQ: circuit with quantum evolution in the chaotic regime; 
q-DJFQ: circuit used as a qubit for measurement of the LE. 



Therefore, we have shown that the DJFQ at high en- 
ergies shows a decay of the LE, similar to other chaotic 



systems [id 1 1 lL Il2l Il3l Il4j . This behavior could be 



experimentally observable if the time scale for the Lya- 
punov decay, TL ya p is much smaller than the decoher- 
ence time Tdecoh due to external sources. We estimate 
r Lyap = h/(r)Ej\) - 0.04-0. 2ns, using A » 0.18 from the 
hard chaos regime and experimental parameters 0, HI- 
In [7J] a value of Td CC oh ~ 20ns was obtained for the 
lowest energy states. At high energies (E ~ 3Ej), we 
find that the spectrum is one order of magnitude more 
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dense than at low energies, therefore Tdccoh should be 
at least one order of magnitude smaller, Tdecoh ~ l ns - 
These roughly estimated values of TL yap ~ 0.1ns and 
T dccoh ~ Ins leave some room for observing the FGR 
and Lyapunov regimes of the LE. However, in order to 
realize an experiment, two issues have to be solved: (i) 
preparation of the initial state and (ii) measurement of 
the fidelity F(t). (i) Preparation of the initial state. In 
order to observe the Lyapunov regime the system has to 
be started from a wave packet narrowly localized in both 
coordinate (phase) and momentum (charge) [l(J 11 1, as 
in Eq.(5). Here we suggest the procedure shown schemat- 
ically in Fig. 4(a). To localize the momentum (charge): 
each of the junctions is connected in parallel to a volt- 
age source V s , which builds up a charge in each junction. 
(For E = 3Ej a voltage of V s « 0.1 mV is needed). To 
localize the coordinate (phase): each of the nodes of the 
DJFQ is connected to a large superconducting reservoir 
through a DC-SQUID as in :23j . Then the system is pre- 
pared staying with V s ^ and "^squid = for a long 
time (t < 0), and at t = the voltage sources are set to 
V s = and the flux in the SQUIDs to $squid = $o/2. 
(ii) Measurement of the fidelity. One has to be able to 
measure the overlap of Eq.(l). One protocol proposed 



23 |26J (see also [L5J) 



originally in [24J and applied in 
consists in coupling the chaotic system under consider- 
ation with a qubit that acts both as a perturbing and 
a measuring device. The approximate hamiltonian for 
this case is H = H <g> |0)(0| + H e ® |1)(1|. H , H e are 
the hamiltonians of the unperturbed and perturbed sys- 
tem, respectively, and |0), |1) are the two basis states of 
the qubit such that when the qubit is in the |1) state it 
induces a perturbation e in the chaotic system. A Ram- 
sey type of experiment is performed: after a time t, ir/2 
pulses are applied to the qubit for two evolutions from 
initial states of the system |\&o) ® (|0) + a\l))/y/2 with 
a = 1 and a = i respectively, from which the fidelity 
amplitude f(t) can be obtained, see [3, 13, 2(| for de- 
tails. A possible implementation of this idea is shown 
in Fig. 4(b). A second DJFQ could be used operating as 



a qubit (called q-DJFQ) for measurement of the LE in 
a DJFQ evolving in the quantum chaotic regime (called 
ch-DJFQ). The q-DJFQ is coupled inductively to the ch- 
DJFQ, such that when the q-DJFQ is in the |0) (|1)) 
state a clockwise (counter-clockwise) current flowing in 
it induces a positive (negative) perturbation in the mag- 
netic flux threaded by the ch-DJFQ. The q-DJFQ should 
have a smaller area than the ch-DJFQ (such that the 
flux induced in the q-DJFQ by the currents in the ch- 
DJFQ is small). For a better observation of the LE the 
ch-DJFQ should be built in a more semiclassical regime 
(r/ ~ 0.1, for example), while the q-DJFQ is built in a 
more quantum regime (j) ~ 0.5, for example). Further- 
more, making measurements with the measuring q-DJFQ 
placed at different distances from the ch-DJFQ could al- 
low to obtain the LE for different perturbation intensi- 



ties. For example, if when placing the q-DJFQ closer to 
the ch-DJFQ the decay of the LE becomes independent 
of the distance, that would indicate that the Lyapunov 
regime was reached. The state preparation and mea- 
surement procedure described here requires coupling to 
several external objects, and thus the resulting decoher- 
ence rates will inevitably increase. However, it has been 
shown in general grounds that TL yap < Td cco h [13]. This 
implies that at least the Lyapunov regime of the LE will 
be observable, and for this reason it is already interest- 
ing to perform the experiment. Moreover, one advantage 
of Josephson nanocircuits is that they can be fabricated 
with well-controlled parameters allowing to study the LE 
for different cases of the effective H (h c g = rf). 
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